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FOREWORD 


This  report  is  presented  in  three  parts;  the  first  part  was  published  as  an  interim  report  for  the 
contract  in  August  98,  and  it  consists  of  a  number  of  theoretical  studies  on  minimization  of  the  g- 
sensitivity  of  strip  resonators  with  various  geometries  and  mounting  structures.  The  second  part  ■ 
describes  an  extension  of  the  theoretical  work,  concentrating  on  some  of  the  blank  geometries 
which  showed  the  highest  probability  of  success,  together  with  practical  results  on  devices.  The 
third  part  describes  the  low-g  sensitivity  SOIL  and  DDIL  OCXOs  development  and  the 
corresponding  test  data. 


VI 


Part  1:  Report  on  The  Design  Optimization  of  Singly  and  Doubly  Rotated  Crystal 
Resonators  Which  Minimize  The  Sensitivity  of  The  Acceleration  Sensitivity  to 
Manufacturing  Imperfections. 


Introduction 

Manufacturing  precision  crystal  resonators  with  low  acceleration  sensitivity  continues  to  be  a 
major  challenge  in  the  production  of  ultra  stable  frequency  sources  for  many  applications,  such 
as  GPS.  This  problem  has  received  quite  a  bit  of  academic  attention  in  recent  years,  and  there 
exist  established  theoretical  results  which  point  to  good  crystal  design  principles  for  low  g- 
sensitivity.  In  production  however,  there  generally  seems  to  be  a  lower  bound  of  a  few  parts  in 
10'''°/g  on  what  is  reproducible  with  any  reasonable  yield,  while  a  small  number  of  crystals  are 
produced  with  very  low  g-sensitivities  (<  1X10*''°/g),  suggesting  that  very  good  resonators  are  at 
least  possible  to  manufacture.  The  fundamental  problem  faced  in  producing  real  resonators  is 
that  the  acceleration  sensitivity  can  be  very  sensitive  to  manufacturing  tolerances  for  a  particular 
design.  The  present  study  seeks  to  quantify  this  phenomenon  by  using  a  finite  element  model  to 
obtain  an  accurate  description  of  the  stresses  and  deformations  in  a  crystal  plate  that  has  been 
mounted  in  various  configurations  and  then  combine  this  solution  with  an  analytical  solution  of 
the  vibration  mode  shape  for  a  harmonic  of  thickness  shear  in  the  general  perturbation 
procedure  developed  by  H.F.  Tiersten  [1,2]  to  obtain  the  frequency  shift.  To  this  end,  a  very 
powerful  software  tool  has  been  created  which  utilizes  previously  developed  finite  element 
software  integrated  together  with  a  custom  developed  numerical  perturbation  integral  solver  to 
compute  the  frequency  shift  under  any  applied  static  bias.  This  finite  element  procedure  has 
been  integrated  into  a  general  nonlinear  optimization  program  to  automate  the  search  for  good 
designs. 

To  compute  the  frequency  shift,  a  general  purpose  integration  method  has  been  developed 
which  calculates  the  first  perturbation  of  the  Eigen-value  from  a  finite  element  solution  of  the 
static  biasing  state  and  an  analytical  solution  of  the  mode  shape  for  a  trapped  energy  resonator 
operating  in  an  arbitrary  overtone  of  thickness  shear.  The  numerical  technique  that  has  been 
employed  generalizes  previously  developed  methods  by  L.D.Clayton  and  E.P.  Eernisse  [5]  with 
many  enhancements  to  improve  accuracy.  Similar  techniques  have  been  applied  to  problems  in 
SAW  acceleration  sensitivity  by  Sinha  and  Locke  [7],  as  well  as  by  the  first  author  [8].  Similar 
finite  element  methods  have  been  developed  for  BAW  devices  by  the  authors  [6]  for  the 
analysis  of  mounting  stresses  in  crystal  resonators.  The  mode  shape  that  has  been  developed 
is  completely  general  accommodating  either  singly  or  doubly  rotated  cuts  of  quartz  through  the 
well  know  planar  transformation  of  the  mode  shape  axes.  The  integration  procedure  readily 
accommodates  this  planar  transformation,  as  well  as  the  existence  of  non-rectangular 
geometries  caused  by  circular  or  elliptical  electrodes  and/or  contouring  by  using  a  very  general 
quadrature  procedure.  Using  the  developed  software  tools,  detailed  sensitivity  studies  are 
performed  on  rectangular,  circular,  elliptical,  and  more  general  shaped  resonator  geometries 
with  various  mounting  schemes  for  the  purpose  of  ranking  the  designs  with  respect  to  their 
sensitivity  to  manufacturing  imperfections. 
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Sensitivity  of  r  To  Misalignment  of  the  Electrode  Plating 

The  sensitivity  of  the  acceleration  sensitivity  to  misalignment  of  the  electrode  is  obtained  by 
following  Tiersten’s  analysis  [3,4],  Shown  in  figure  (2)  is  the  geometry  of  a  rectangular  crystal 
resonator  with  a  single  pair  of  identical  electrodes,  which  in  the  unperturbed  case,  are  located  in 
the  center  of  the  plate.  When  the  resonator  is  manufactured  it  is  assumed  to  have  a  small  error 
in  the  electrode  placement,  as  shown.  The  resulting  electrode  will  be  displaced  from  the  center 
by  a  small  distance  R,  at  an  angle  a.  In  general,  the  acceleration  sensitivity  vector  for  the 
resonator  is  denoted  by 


(1) 


The  sensitivity  of  r  to  electrode  misalignment  is  obtained  by  expanding  this  quantity  to  first 
order  in  the  mode  offset  distances 

'  '  (2) 


'U) 


+  • 


dd. 


-d^,A  =  l,3 


where  d,  =Rcosa  and  d^  =i?sinor  are  the  mode  offset  distances  in  the  x,  and  X3  directions, 
respectively.  For  a  perfectly  symmetric  mounting  structure. 


r<'>  =  0  (3) 

which  gives 

^  r<'>  =  cy  V,  (4) 

where  Ca^  are  the  so-called  mispositioning  coefficients.  The  maximum  sensitivity  to  mode 
mispositioning  is  defined  as 

(5) 


where 

Tmax  =  inax{|r|  ,a  =  0°...360°}  (6) 

with 


|r|  =  7r^')2^p(2)2  ^  p(3)2 

representing  the  magnitude  of  the  acceleration  sensitivity  vector.  Cmax  represents  the  maximum 
sensitivity  of  the  acceleration  sensitivity  to  misalignment  of  the  electrode  structure  with  respect 
to  the  support  system. 
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Minimization  of  the  Sensitivity  of  r  to  Manufacturing  Imperfections  as  a  Nonlinear 
Optimization  Problem 

For  a  given  support  structure,  the  mode  mispositioning  coefficient  is  a  function  of  geometric  and 
material  parameters  defining  the  structural  system  of  resonator  and  mount.  In  general, 


(8) 

where  a,-  represent  N  generalized  variable  parameters  defining  the  geometry,  crystal  cut,  blank 
orientation,  etc.  Define  the  general  nonlinear  optimization  problem: 

Minimize 

Subject  to:  af  <  a.  <  (9) 

where  af  and  represent  lower  and  upper  bounds  placed  on  each  parameter. 

Figure  (1)  shows  a  functional  layout  of  the  non-linear  optimization  program  developed  to  solve 
equation  (9).  The  present  analysis  makes  use  of  the  IMSL  routine  BCONF,  which  minimizes  a 
function  of  N  variables  subject  to  bounds  on  the  variables  using  a  quasi-Newton  method  and  a 
finite  difference  gradient. 

Discussion  of  Results 

Using  the  developed  software,  numerical  studies  have  been  carried  out  on  resonators  of 
different  geometries  and  crystal  cuts  for  the  purpose  of  comparing  the  various  designs.  In  all 
cases  the  plate  thickness  has  been  held  constant  and  a  fundamental  thickness  harmonic  has 
been  considered.  Also,  the  basic  ranges  of  major  dimensions  have  been  held  constant  for 
different  blank  geometries  and  mounting  schemes.  This  allows  the  different  designs  to  be 
compared  most  accurately.  The  basic  dimensions  used  are  summarized  in  figure  (2)  and  figure 
(3)  shows  the  sense  of  rotation  for  the  blank  rotation  used  in  ail  the  studies.  Figures  (4. a)  -  (4.j) 
show  the  results  of  a  basic  study  of  the  effect  of  a  0.1  mm  mode  offset  on  the  total  r  of  a 
rectangular  SC  Cut  resonator  rigidly  mounted  along  all  four  edges.  These  plots  show  the 
magnitude  of  r  in  units  of  10'’°/g  vs.  the  angle  of  mode  offset,  a,  for  various  length  to  width 
ratios.  This  study  was  done  for  LTW  ratios  of  0.5,  0.6,  0.7,  0.8,  0.9,  1.0,  1.11,  1.25,  1.43,  1.67, 
and  2.0.  Only  the  curves  representing  L/W  =  0.5,  1.0,  and  2.0  are  labeled  so  as  to  not  clutter  the 
graphs.  Each  of  the  graphs  in  figure  (4)  represent  different  blank  rotations  with  respect  to  the 
crystallographic  axes,  taken  in  twenty  degree  increments.  From  these  curves  it  is  observed  that 
the  sensitivity  of  the  resonator  to  mode  offset  varies  widely  with  respect  to  blank  orientation  and 
planar  aspect  ratio  with  roughly  an  order  of  magnitude  difference  between  the  best  and  worst 
values.  Figure  (4.k)  summarizes  the  results  by  plotting  Cmax  as  a  function  of  planar  aspect  ratio 
for  different  blank  orientations.  The  horizontal  scale  in  this  graph  is  in  dimensionless  aspect  ratio 
with  the  vertical  scale  representing  units  of  10"®  per  meter  per  1  g  acceleration.  By  tracing  the 
minimum  Cmax  over  all  aspect  ratios,  it  can  be  seen  that  for  this  particular  mounting  condition, 
the  1:1  aspect  ratio  is  the  worst  case  with  the  best  cases  realized  by  elongated  plates  (1:2  or 
2:1  planar  aspect  ratios).  Figure  (6)  shows  a  refined  study  done  in  the  vicinity  of  y=136®  to  156° 
which  shows  an  approximate  minimum  at  7=128°.  This  result  was  verified  later  with  the 
optimizer.  Figure  (5)  shows  similar  results  for  larger  aspect  ratios.  It  is  observed  here  that  the 
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sensitivity  tends  to  level  off  at  higher  aspect  ratios,  suggesting  that  for  greatly  elongated  plates, 
the  Cmax  is  not  sensitive  to  errors  in  actual  blank  dimensions.  Figure  (7)  shows  the  dependency 
of  Cmax  on  y  for  the  range  of  aspect  ratios  from  1:2  to  2:1  for  an  SC  Cut  rectangular  plate  rigidly 
mounted  along  all  four  edges.  Here  it  is  seen  again  that  the  optimum  blank  orientation  occurs  at 
38°  +/-  90°.  Figure  (8)  shows  the  results  of  the  same  study  carried  out  on  a  plate  with  rigid 
mounts  along  only  two  edges  (xi  =  +/-  L).  Here  it  is  seen  that  the  best  result  is  obtained  for  an 
aspect  ratio  of  1 :2  at  roughly  the  same  angle.  This  represents  the  case  of  a  rigid  mount  along 
the  long  dimension  of  the  plate.  The  case  of  a  rigid  mount  along  the  short  dimension  is 
represented  by  aspect  ratios  greater  than  1.0.  For  this  case  the  optimum  occurs  at  a  slightly 
higher  value  of  y  with  about  twice  the  sensitivity  of  that  obtained  for  the  four  edge  mount.  Figure 

(9)  shows  the  results  of  a  similar  study  carried  out  for  an  FC  Cut  plate.  By  inspection  it  is  seen 
that  the  FC  Cut  behaves  essentially  the  same  as  the  SC  with  only  slight  differences.  Figures 

(10)  -  (12)  represent  a  study  carried  out  for  an  AT  Cut  plate  and  figure  (13)  shows  results  for  the 
BT  Cut.  It  is  observed  that  the  theoretical  behavior  of  the  singly  rotated  cuts  is  very  different 
from  that  of  the  SC  or  FC.  The  optimum  values  obtained  for  AT  and  BT  cuts  are  found  at  y  =  0° 
with  roughly  twice  the  sensitivity  of  the  four  edge  mounted  SC  Cut.  Shown  in  figure  (14)  is  the 
geometry  of  a  rectangular  SC  Cut  resonator  mounted  along  two  edges  with  thin  slotted  ribbons 
with  overhang  of  the  crystal  beyond  the  support  by  a  distance  of  U.  Figures  (14.a)  -  (14.e)  show 
the  results  of  a  study  carried  out  on  this  plate  for  various  values  of  U-  Considering  the  case  of 
mounting  along  the  short  edge  (LA/V  >  1),  it  is  observed  here  that  the  minimum  Cmax  is  found  at 
roughly  y=  120°  for  all  cases  with  the  actual  minimum  value  changing  with  overhang  distance. 
Using  the  optimization  program,  optimal  values  have  been  found  for  this  case  of  L  =  5.0  mm,  W 
=  2.9  mm  (L7W  =  1.72),  y  =  129°,  and  U  =  3.5  mm  (U/Ltotai  =  0,206).  Figures  (15.a)  and  (15.b) 
show  the  centerline  deflection  of  the  blank  along  the  two  planar  axes  under  the  action  of  normal 
acceleration  for  different  U  values.  Interesting  to  note  here  is  that  the  optimum  value  does  not 
correspond  to  the  minimum  deflection  case.  Figure  (15.c)  shows  the  effect  of  overhang  distance 
on  Cmax-  The  dotted  line  in  this  figure  represents  the  same  plate  rigidly  mounted  along  the  two 
edges  with  the  same  inner  dimensions.  Figures  (16)  -  (18)  show  results  obtained  for  an  elliptical 
SC  Cut  plate  using,  respectively,  two  point,  four  point,  and  fully  supported  mounting  schemes. 
Interesting  to  note  is  that  in  the  fully  supported  case,  there  is  very  little  difference  between  the 
elliptical  plate  and  the  rectangular  plate  as  shown  in  figure  (7).  Also  worth  noting  is  that  in  the 
case  of  the  four  point  mount  as  shown  in  figure  (17),  the  optimal  values  are  obtained  at  y  =  0°, 
+/-  90°.  Figure  (19)  shows  the  results  of  a  study  done  on  the  effect  of  geometric  scale  on  Cmax 
for  an  SC  Cut  rectangular  plate  rigidly  supported  along  all  four  edges.  In  this  figure,  the 
magnification  of  the  Cmax  value  is  plotted  as  a  function  of  geometric  scale  factor,  K,  for  various 
cases.  From  these  results,  it  is  observed  that  Cmax  is  relatively  invariant  to  scale  when  the 
thickness  of  the  plate  is  scaled  with  the  lateral  dimensions,  but  is  very  sensitive  to  lateral  scaling 
at  constant  thickness.  In  general,  Cmax  will  scale  with  lateral  dimensions  at  constant  thickness 
(constant  frequency)  as  the  square  of  K  when  the  electrode  size  is  held  constant.  Figure  (20) 
shows  the  configuration  for  a  general  shape  blank  in  which  a  parametric  cubic  spline  is  used  to 
define  the  outer  edge  of  the  plate.  Shown  in  figures  (20.a)  and  (20.b)  are  the  possible  shapes 
realizable  by  the  curves  when  the  shape  factors  ki  and  ka  are  varied  through  values  of  1.0  to  3.0 
for  convex  and  concave  profiles.  Figures  (21. a)  -  (21  .j)  show  the  results  of  a  study  done  for 
convex  and  concave  geometries  using  various  ki  and  ka  values  for  a  four  point  mounted  SC  Cut 
resonator.  It  is  observed  here  that  while  the  general  shape  of  the  plotted  curves  changes  with 
the  geometry,  the  minimum  Cmax  obtainable  does  not  vary  appreciably.  This  fact  is  highlighted  in 
the  plots  of  figure  (23)  which  show  the  effect  of  varying  the  shape  factors  ki  and  ka  on  the 
optimum  Cmax-  The  total  variation  in  Cmax  shown  in  these  plots  is  roughly  0.1X10"®  /m/g,  which  is 
quite  small  when  compared  to  the  effect  of  aspect  ratio  and  blank  orientation. 


4 


Part  One  Conclusions 


Shown  in  figure  (23)  is  a  table  of  results  obtained  from  the  optimization  program  for  the  various 
geometries  and  crystal  cuts  considered.  The  very  best  results  are  obtained  from  the  fully 
supported  plate  with  crystal  cut  in  the  vicinity  of  the  SC,  with  Cmax»  0.25.10'®.  This  value  is 
achieved  at  a  specific  aspect  ratio  and  blank  orientation,  and  is  generally  not  changed 
dramatically  by  using  different  shapes  other  than  a  rectangle.  The  (f)  angle  has  a  small  effect  on 
the  minimum  Cmax,  although  this  is  not  a  dominating  factor.  In  the  case  of  the  rectangular  plate 
rigidly  mounted  along  two  edges  only,  it  is  observed  that  the  minimum  value  of  Cmax=0.56.10'® 
for  the  SC  Cut  is  roughly  twice  that  obtained  for  the  full  four  edge  mount.  For  the  AT  Cut,  an 
optimum  value  of  Cmax=0.50.10'®  is  obtained  for  the  four  edge  mount  and  a  corresponding  value 
of  Cn,ax=0.80.10'®  is  obtained  for  the  two  edge  mount.  It  appears  that  the  compensating  effects 
of  optimal  blank  aspect  ratio  and  orientation  are  fully  developed  with  a  continuous  mount  around 
the  perimeter.  For  a  two  edge  mount,  the  effect  of  overhangs  has  been  found  to  be  a  factor  at  a 
given  blank  orientation  as  illustrated  by  figure  (15.c).  In  the  example  considered  here,  the 
overhang  tends  to  improve  the  sensitivity  of  the  structure  to  manufacturing  imperfections  over 
the  corresponding  rigid  mounted  plate.  The  effects  of  scaling  the  geometry  as  illustrated  by 
figure  (19)  are  also  significant.  These  results  suggest  that  an  advantage  can  be  gained  by 
shrinking  the  planar  geometry  of  a  crystal  at  a  given  frequency.  For  a  four  point  mounted  crystal, 
there  is  a  small  effect  of  overall  shape  on  the  minimum  Cmax  obtainable.  The  results  of  the 
shape  optimization  suggest  that  there  is  a  small  advantage  to  be  gained  by  effectively  trimming 
off  the  portions  of  an  elliptical  blank  between  the  supports,  resulting  in  somewhat  of  a  diamond 
shape.  In  general,  the  most  important  parameters  for  reducing  Cmax  are  the  planar  aspect  ratio 
and  the  blank  orientation,  and  not  the  actual  shape  of  the  resonator. 
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Figure  1.)  Functional  layout  of  non-linear  optimization  software  package. 


A  s’3 


^ _ 

_ 

2a 

R  =  Mode  offset  distance 
a  =  Mode  offset  angle 


2L 


2a  =  4.0  mm 
2b  =  4.0  mm 
2L  =  5.0  mm  -  10.0mm 
2w  =  5.0  mm  -  10.0mm 
2h  =  0.1644  mm 
2h’  =  50  nm  An 

Figure  2.)  Geometry  and  dimensions  for  mode 
offset  study. 


y  =  Blank  rotation  angle 
2h  =  Blank  thickness 
2h’  =  Electrode  thickness 


Figure  3.)  Orientation  of  crystal  blank 
rotation  for  mode  offset  study. 
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Figure  4.)  Sensitivity  study  for  SC  Cut 


Figure  4.)  Sensitivity  study  for  SC  Cut 


Figure  4.)  Sensitivity  study  for  SC  Cut 
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Figure  8.)C„3^  vs.  y  for  SC  Cut  at  various  LAV  values  for  rigid  mount 

along  two  edges  (x-x). 
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Figure  9.)  Sensitivity  study  for  FC  Cut 


Figure  9.)  Sensitivity  study  for  FC  Cut 


Figure  10.)  Sensitivity  study  for  AT  Cut 
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Figure  10.)  Sensitivity  study  for  AT  Cut 


Figure  13.)  Sensitivity  study  for  BT  Cut 


Figure  13.)  Sensitivity  study  for  BT  Cut 


Figure  13.)  Sensitivity  study  for  BT  Cut 
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Figure  14.)  Mode  mispositioning  analysis  for  slotted  ribbon  mounted  rectangular  SC  Cut 
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Figure  14.)  Mode  mispositioning  analysis  for  slotted  ribbon  mounted  rectangular  SC  Cut 


d.)  Lq  =  1.00  mm 
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Figure  14.)  Mode  mispositioning  analysis  for  slotted  ribbon  mounted  rectangular  SC  Cut 

plate. 
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Figure  18.)  Mode  mispositioning  analysis  for  continuously  mounted  elliptical  SC  Cut  plate 
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Figure  19.)  Effect  of  geometric  scale  on  the 
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ositioning  coefficient  for  a  rectangular  SC  Cut  plate. 
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Xi(t)  =  (2t3-3t2+l)Bj  +  (-2t3-3t2)L,  +  (t3-2t2+t)  kjL, 
X3(t)  =  (2t3-3t2+l)L3  +  (-2tH3t2)  B3  -  k3L3 
for  convex  shape 

x,(t)  =  (2t3-3t2+l)Bi  +  (-2t3-3t2)Li  +  (t3-t2)  k3L3 

X3(t)  =  (2t3-3t2+l)L3  +  (-2t3+3t2)  B3  -  (t3-2t^t)  kjLj 

for  concave  shape 


Cmax  X  10  /m/g  Cmax  X  10  /m/g 


Blank  Orientation  (deg.) 
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Figure  21.)  vs.  blank  orientation  for  convex  and  concave  blank  geometries. 
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OraxXIO  /m/g  CiraxXIO  /m/g 
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Figure  21.)  vs.  blank  orientation  for  convex  and  concave  blank  geometries. 


33 


e.)  Convex  Shape ,  kj  =  kj  =  3.0 


CitBxXIO  /m/g  CmaxXIO  /m/g 


g.)  Concave  Shape ,  kj  =  kj  =  0.25 
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Blank  Orientation  (deg.) 
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Blank  Orientation  (deg.) 

Figure  21.)  Cn,aj  vs.  blank  orientation  for  convex  and  concave  blank  geometries. 
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i.)  Concave  Shape,  kj  =  kj  =  0.75 


Cmax  at  Optimal  Parameters  Cmax  at  Optimal  Parameters 


a.)  Convex  Profile 


0.0  0.2  0.4  0.6  0.8  1.0  1.2 


ki  =  ka 


Figure  22.)  Minimum  C„ax  shape  factor  for  convex  and  concave  blank  geometries. 
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Figure  23.)  Summary  of  results  from  non-linear  optimization  study. 
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Part  2:  Report  on  The  Theoretical  and  Experimental  Results  for  Optimization  of 
the  Acceleration  Sensitivity  of  Rectangular  Crystal  Resonators 


Introduction 

Theoretical  and  experimental  results  for  the  acceleration  sensitivity  of  rectangular  AT  and  SC 
Cut  crystal  resonators  are  presented.  In  a  previous  paper  [1],  a  detailed  theoretical  analysis  was 
presented  on  the  g-sensitivity  of  crystal  resonators  with  non-circular  geometry.  This  work  was 
also  described  in  detail  in  the  interim  report  [2],  which  is  also  presented  in  more  detail  as  the 
first  part  of  this  report. 

The  present  study  focuses  on  the  realization  of  various  designs  taken  from  the  first  analysis. 
Rectangular  AT  Cut  crystal  resonators  have  been  fabricated  and  mounted  in  a  rigid  two-ended 
configuration  and  a  two-ended  with  overhang  configuration.  The  effect  of  support  spacing  on  the 
g-sensitivity  for  these  crystals  is  discussed  and  the  results  are  compared  with  conventional  rigid 
mounting  (no  overhang).  Rectangular  SC  Cut  crystals  have  also  been  fabricated  and  the  effect 
of  blank  orientation  on  the  g-sensitivity  for  these  crystals  is  presented.  A  brief  discussion  of  the 
g-sensitivity  measurement  technique  is  also  presented. 

Manufacturing  precision  crystal  resonators  with  low  acceleration  sensitivity  continues  to  be  a 
major  challenge  in  the  production  of  ultra  stable  frequency  sources  for  many  applications,  such 
as  GPS.  The  fundamental  problem  faced  in  producing  real  resonators  is  that  the  acceieration 
sensitivity  can  be  very  sensitive  to  manufacturing  tolerances  for  a  particular  design. 
Misplacement  of  the  electrode  plating,  misalignment  of  the  crystal  blank  in  the  supports,  and 
inaccuracies  in  the  amount  and  distribution  cf  adhesive  all  lead  to  appreciable  degradation  in 
the  g-sensitivity  of  manufactured  resonators.  In  a  previous  paper  [1],  theoretical  results  were 
presented  for  the  analysis  of  the  acceleration  sensitivity  of  singly  and  doubly  rotated  crystal 
resonators  with  non-circular  geometry.  The  results  of  this  study  showed  that  a  substantial 
improvement  could  be  obtained  by  using  crystal  blanks  with  rectangular  geometry,  and  more 
recent  results  have  shown  that  there  are  definite  optimal  blank  dimensions  for  particular 
mounting  schemes.  In  general,  it  has  been  found  that  the  major  effects  are  the  planar  aspect 
ratio  and  orientation  of  the  plate.  The  analytical  studies  have  also  revealed  that  blank  size 
reduction  at  a  given  frequency  and  the  use  of  interior  supports  which  result  in  “overhangs”  on 
the  mounted  crystal  may  also  play  a  role  in  improving  the  g-sensitivity  of  crystal  resonators. 
Based  on  the  developed  design  criteria,  AT  and  SC  cut  rectangular  resonators  have  been 
fabricated  using  a  variety  of  mounting  strategies  and  blank  orientations.  Experimental  results  for 
these  crystals  are  presented  along  with  the  corresponding  theoretical  predictions. 

Measurement  Technique 

The  setup  shown  in  figure  (1)  is  used  for  measurement  of  discrete  resonators.  The  method  is  a 
variant  of  the  technique  reported  by  Filler  [3]  and  others,  in  which  a  crystal  oscillator  is 
mechanically  mounted  on  a  vibration  table,  and  a  low  frequency  sinusoidal  vibration  field  is 
applied.  The  effect  is  simple  frequency  modulation  with  a  very  low  modulation  index,  and  the 
relationship  between  the  sideband  levels  and  the  g-sensitivity  is  shown  to  be: 


acoo 
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where  L  is  the  sideband  level  in  dBc,  y  is  the  g-sensitivity,  a  is  the  g  level,  ax,  is  the  vibration 
frequency,  and  axo  is  the  oscillator  frequency. 

In  the  case  of  discrete  crystal  measurement,  the  crystal  is  physically  mounted  on  the  vibration 
table  with  a  thin  coaxial  cable  connecting  it  to  a  180°  hybrid  transformer.  On  the  balancing  port 
of  the  hybrid  a  matched  cable  length  with  capacitive  trimmer  is  used  to  cancel  out  the  Co  of  the 
resonator  and  the  connecting  cable.  This  is  a  convenient  network  as  it  gives  a  frequency 
response  which  is  symmetrical  about  the  series  resonance.  This  is  a  passive  network,  and  the 
signal  source  is  an  external  synthesizer  (in  our  case  we  used  an  HP3335)  set  to  the  crystal's 
series  resonance.  The  relationship  between  the  sideband  levels  and  the  g-sensitivity  is  now 
defined  by  the  equation: 

where  0x3  is  the  3  dB  bandwidth  of  the  crystal  response. 

For  a  vibration  environment  for  this  test,  a  frequency  of  70  Hz  has  been  found  to  be  convenient, 
with  a  level  of  lOg  peak  acceleration.  This  avoids  interaction  with  60  Hz  or  50  Hz  mains 
frequencies  and  gives  a  workable  displacement  level.  Lower  frequencies  can  be  used  to  good 
effect  to  increase  the  vibration  induced  sideband  levels,  but  this  can  result  in  rather  high 
mechanical  displacements.  A  fixture  has  been  made  to  enable  the  devices  to  be  conveniently 
measured  in  three  mutually  orthogonal  axes. 

The  HP4396  Network  /  Spectrum  Analyzer  is  particularly  well  suited  to  this  measurement  as  it  is 
one  of  a  genera  of  instruments  which  performs  its  IF  filtering  by  an  FFT  technique.  This  enables 
it  to  operate  with  relatively  short  sweep  times  even  with  resolution  bandwidths  as  low  as  1  Hz. 
The  measurement  sequence  is  set  up  in  the  instrument's  internal  programming  language,  and 
starts  with  the  analyzer  in  'Network'  mode  and  the  synthesizer  switched  off.  Successive  sweeps 
are  used  to  determine  the  center  frequency  and  bandwidth  of  the  resonator.  Then  the 
instrument  is  switched  to  'Spectrum'  mode,  and  the  synthesizer  is  set  to  the  center  frequency  of 
the  response.  For  a  70  Hz  vibration  frequency,  a  span  of  300  Hz  is  set  up  on  the  spectrum 
analyzer,  centered  on  the  synthesizer  frequency,  with  a  resolution  bandwidth  of  1  Hz.  The 
instrument  markers  are  used  to  determine  the  sideband  levels  relative  to  the  carrier,  and  this 
information  is  used  to  calculate  the  g-sensitivity  directly. 

For  some  devices,  motional  parameters  were  also  characterized  using  a  Transat  CNA  test 
station.  This  system  uses  a  pi-network  as  a  measurement  port  and  performs  impedance 
measurements  in  a  pre-defined  sequence  around  the  crystal  resonance.  These  data  are 
manipulated  using  the  method  described  in  the  IEC444-5  measurement  standard  to  derive  the 
Co,  Cl,  Ri  and  Q  of  each  resonator. 

Discussion  of  Results 

The  first  study  consists  of  a  rectangular  AT  Cut  resonator  mounted  into  a  Vectron  S-Type 
ceramic  package,  as  shown  in  figure  (2).  This  package  has  shelves  at  two  opposite  ends,  which 
allow  for  double-ended  mounting.  These  crystals  were  fabricated  with  a  rectangular  electrode 
that  is  centered  with  respect  to  the  crystal  blank  and  mount.  These  crystals  were  directly 
mounted  to  the  ceramic  with  manually  dispensed  conductive  epoxy  adhesive.  Shown  in  figure 
(3)  are  the  results  of  an  experiment  performed  on  this  configuration  for  different  blank 
frequencies.  In  this  experiment  a  single  ended,  or  cantilever,  mount  was  used.  As  expected,  the 
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component  of  r  normal  to  the  crystal  plate  dominates  the  g-sensitivity  of  these  devices  vAth  a 
nominal  value  of  about  3x10'®/g.  Interesting  to  note  from  this  plot  is  how  the  spread  in  the  data 
tends  to  increase  (particularly  above  30  MHz),  while  the  average  value  remains  relatively 
constant  at  3x10'^/g.  As  the  frequency  increases,  the  plate  thickness  decreases,  and  it  appears 
that  the  sensitivity  of  r  to  manufacturing  imperfections  likewise  increases.  The  likely  source  of 
manufacturing  imperfections  in  these  devices  is  the  inaccuracy  of  adhesive  distribution  inherent 
in  the  mounting  process.  At  30.0  MHz  and  R  =  2.5  %,  an  AT  Cut  crystal  is  roughly  0.0021  in 
thick,  and  at  40.0  MHz,  this  value  drops  to  0.0016  in.  Considering  a  typical  adhesive  thickness 
of  0.001  in  -  0.002  in,  it  appears  that  as  the  crystal  thickness  approaches  the  adhesive 
thickness,  the  mispositioning  sensitivity  effect  becomes  substantial  and  the  spread  in  the  results 
increases. 

Figures  (4)  -  (7)  show  the  results  of  a  study  carried  out  on  a  16.384  MHz  rectangular  AT  Cut 
crystal  mounted  in  the  same  configuration.  A  series  of  crystals  were  built  in  both  the  X’  and  Z’ 
elongated  configurations,  and  mounted  with  both  single  and  double-ended  methods.  The  results 
obtained  show  very  different  behavior  for  the  four  different  configurations.  The  best  results  were 
obtained  from  the  Z’-elongated  crystals  with  double  ended  mounting,  resulting  in  nominal  r 
values  at  5X10‘^°/g  and  a  spread  of  roughly  2X10'^°/g.  The  worst  cases  were  found  with  the  X’- 
elongated  crystals.  The  X’-elongated  devices  with  double-ended  mount  showed  nominal  values 
of  2X10'®/g  -  3X10'®/g  with  a  huge  spread  in  the  data.  One  plausible  explanation  for  these 
differences  is  the  difference  in  energy  trapping  found  in  these  two  orientations.  The  AT  Cut  traps 
better  along  the  Z’  direction  than  it  does  along  the  X’  direction.  This  means  that  for  the  Z’- 
elongated  plate  the  mode  shape  is  better  centered  and  center-localized  which  leads  to  good 
symmetry  and  less  interaction  of  acoustic  energy  with  the  mounts. 

Figure  (8)  shows  the  results  of  a  study  carried  out  on  very  small  rectangular  crystals.  A  set  of 
0.140  in  X  0.060  in  crystals  used  in  the  Vectron  K-Type  oscillator  were  mounted  using  four 
different  schemes.  These  crystals  showed  relatively  low  numbers  in  the  1X10'®/g  -  5X10'®/g 
range,  but  did  not  exhibit  any  meaningful  trend.  The  results  here  show  a  relatively  large  spread 
with  no  noticeable  differences  between  different  mounting  methods.  These  crystals  were 
mounted  to  ceramic  using  conductive  epoxy  adhesive.  It  is  likely  here  that  the  distribution  of 
adhesive  plays  a  role  in  the  scattering  of  the  data. 

Figures  (9)  -  (12)  show  the  results  of  a  theoretical  study  done  with  the  previously  developed 
software  [1].  In  this  study,  a  rectangular  crystal  was  considered  as  continuously  mounted  onto  a 
pair  of  supporting  bars  with  epoxy  adhesive,  resulting  in  an  overhang  on  each  side.  The 
calculations  were  performed  for  various  beam  spacings  and  overhang  lengths  for  a  plate  width 
of  0.076  in.  and  a  frequency  of  20  MHz.  As  seen  in  figure  (11),  there  are  optimal  dimensions  for 
which  the  g-sensitivity  of  this  configuration  can  be  made  very  small.  Figure  (12)  shows  the 
results  of  a  similar  study  performed  on  a  0.210  in  X  0.076  in  Z’-elongated,  20MHz  crystal.  In 
figure  12.a  the  actual  r  components  are  shown  separately,  and  it  is  noted  that  the  normal 
component  of  the  g-sensitivity  changes  sign  at  a  beam  spacing  of  roughly  0.125  in.  This  is  a 
very  important  point  because  it  means  that  the  normal  g-sensitivity  can  be  made  arbitrarily  small 
by  optimizing  the  beam  spacing  for  a  particular  crystal. 

To  test  this  hypothesis,  a  series  of  rectangular  AT  Cut  crystals  were  fabricated  and  mounted 
onto  a  lead  frame  structure  as  shown  in  figure  (13).  The  lead  frame  was  photo-chemically 
etched  from  stainless  steel  with  pull-tabs  which  are  removed  after  bonding  to  the  substrate.  This 
structure  provides  a  very  accurate  geometry,  and  the  pull-tab  approach  enables  connections  to 
be  easily  made  to  the  crystal  electrodes.  Once  the  pull-tabs  have  been  removed,  a  crystal  with 
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dual-ended  termination  of  metallization  is  mounted  with  conductive  epoxy  as  shown  in  figure 
13.b.  For  this  purpose,  a  base-plating  mask  was  developed  with  alignment  marks  to  help  center 
the  blank  onto  the  mount. 

Figures  (14)  -  (16)  show  the  results  obtained  from  these  crystals  for  beam  spacing  of  0.120, 
0.130,  0.140,  and  0.150  in.  It  is  observed  here  that  the  trend  in  the  data  is  towards  an  optimal 
value  of  0.140  in  -  0.150  in,  which  follows  the  theoretical  data  in  principle,  with  a  shift  in  the  best 
value  from  the  0.125  in  prediction.  It  should  also  be  noted  from  figure  (15)  that  the  in-plane  r 
components  remain  relatively  constant  at  1-2X1 0‘''°/g.  Figure  (16)  shows  a  comparison  of 
results  obtained  by  driving  fundamental  and  3'**  overtone  modes  in  the  same  devices.  The  third 
overtone  results  reveal  a  slight  minimum  at  the  0.140  in  spacing.  It  was  also  found  in  these 
crystals  that  the  spread  in  results  for  the  3"*  overtone  was  tighter  than  that  of  the  fundamental. 
This  is  likely  caused  by  the  superior  energy  trapping  found  in  the  3"'*  overtone  as  compared  with 
the  fundamental. 

Based  on  these  preliminary  results,  two  larger  groups  of  resonators,  a  total  of  120  units,  were 
carefully  manufactured  to  the  same  design  with  a  0.140”  beam  spacing.  One  group  was  plated 
to  optimize  the  resonator  performance  at  20  MHz  on  the  fundamental  mode  and  the  second 
group  was  optimized  for  third  overtone  operation  at  60  MHz.  Initially  these  devices  were 
measured  with  the  passive  technique,  but  the  absolute  measurement  resolution  with  this 
method  is  in  the  order  of  1.10'''°  for  20  MHz  resonators  with  Q  factors  below  100,000.  Since 
many  of  the  units  are  well  below  this  value  in  at  least  one  axis,  oscillators  were  made  for  use  on 
the  vibration  fixture.  For  this  type  of  resonator  this  improves  the  measurement  resolution  by  an 
order  of  magnitude.  The  results  are  summarized  below: 

20  MHz  Fundamental  Units 


Q  factor 

Ri 

Co 

Ci 

Fi 

Fa 

Fo 

|F| 

73.9 

19.0 

4.2 

6.32 

0.09 

0.11 

0.22 

a 

18.7 

9.1 

0.9 

0.17 

0.05 

0.06 

0.08 

Min 

19.7 

11.9 

3.7 

0.01 

0.01 

0.03 

0.09 

Max 

106.3 

67.1 

7.4 

6.61 

0.20 

0.38 

0.27 

0.46 

60  MHz  3"*  Overtone  Units 


Q  factor 

Ri 

Co 

Ci 

Fi 

F2 

Fo 

|F| 

Mean 

47.5 

102.2 

3.9 

0.59 

0.08 

0.16 

0.17 

0.27 

a 

11.2 

27.1 

0.1 

0.05 

0.06 

0.11 

0.09 

0.10 

Min 

27.6 

62.3 

3.5 

0.46 

0.01 

0.01 

0.02 

0.10 

Max 

67.7 

162.3 

4.1 

0.72 

0.25 

0.36 

0.63 

In  these  results.  Fa  is  the  normal  component.  and  Fa  are  the  in-plane  components,  with  Fa 
along  the  short  axis  and  Fi  along  the  length  of  the  blank.  The  full  distributions  are  tabulated  in 
Appendix  A  and  shown  graphically  in  fig  17  a  to  d.  It  is  clear  from  these  results  that  it  is  practical 
to  obtain  g-sensitivity  results  for  this  class  of  resonator  in  the  low  parts  in  10^°  per  g,  and  that  in 
this  case  the  choice  of  overtone  is  not  a  dominant  parameter. 

Figure  (18)  shows  the  results  of  a  study  done  on  rectangular  SC  Cut  resonators.  In  this  study,  a 
series  of  rectangular  20  MHz  SC  Cut  resonators  with  nominal  dimensions  of  0.220  in  X  0.090  in 
were  fabricated  by  dicing  the  blanks  from  larger  wafers  at  oblique  orientations  with  respect  to 
the  Z’  axis  of  the  crystal.  These  crystals  were  mounted  into  the  ceramic  S-Type  package  as 
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shown  in  figure  (2)  and  the  r  components  measured.  The  results  obtained  from  these  crystals 
show  a  wide  variation  with  the  best  values  obtained  at  0,  22.5,  and  90°.  The  worst  values  seen 
in  these  devices  were  obtained  at  -22.5,  -45.0,  and  -67.5°.  At  the  best  orientations  nominal 
values  of  less  than  1X10'®/g  were  obtained,  but  overall  the  performance  of  these  crystals  was 
found  to  be  poorer  than  the  equivalent  AT  Cut  devices. 

The  final  study  to  discuss  is  shown  in  figures  (19)  -  (21).  In  this  study,  a  series  of  slotted  clips 
were  fabricated  by  photochemical  etching  of  0.003  in  thick  stainless  steel,  as  shown  in  figure 
(19).  Using  these  clips,  a  rectangular  20  MHz  AT  Cut  Z’-elongated  crystal  was  mounted  to  a 
ceramic  substrate  with  small  mounting  blocks  as  shown  in  figure  (20).  Several  devices  were 
fabricated  at  different  clip  spacings,  and  these  results  are  shown  in  figure  (21).  As  with  the  lead 
frame  devices,  the  in-plane  F  components  are  constant  at  roughly  1-2X1 0'’°/g,  while  the  normal 
components  vary  with  support  spacing.  In  general,  the  results  obtained  with  these  devices  were 
not  as  good  as  those  obtained  with  the  lead  frames.  It  had  been  expected  that  the  improved 
symmetry  through  the  thickness  obtained  with  the  slotted  clips  would  improve  the  g-sensitivity 
(particularly  the  in-plane  components)  of  these  devices,  it  is  possible  that  the  high  profile 
obtained  with  the  mounting  blocks  in  combination  with  their  relatively  large  mass  degrades  the 
performance  of  these  devices.  Future  efforts  with  the  slotted  clip  mounts  should  focus  on  lower 
profile  mounts. 

Part  Two  Conclusions 

Very  low  g-sensitive  resonators  can  be  manufactured  using  rectangular  AT  Cut  Z’-elongated 
crystals  with  a  double-ended  mount  that  incorporates  a  small  overhang  of  the  crystal  beyond  the 
mount.  For  a  well  made  device  in  this  configuration,  the  in-plane  components  will  remain 
constant  at  roughly  1-2X10'^°/g  while  optimal  crystal  length  and  beam  spacing  can  adjust  the 
normal  component.  It  appears  that  devices  of  this  type  can  be  reliably  fabricated  with  a  total  □ 
of  2.0-3. 0X1 0‘^°/g  for  fundamental  or  3*^^  overtone  with  a  reasonable  yield.  To  produce  better 
devices,  methods  for  reducing  the  in-plane  components  must  be  addressed.  The  most  obvious 
way  to  achieve  this  is  by  improving  the  symmetry  of  the  mount  in  the  thickness  direction. 
Furthermore,  crystal  attachment  methods  must  be  improved  to  minimize  the  spread  in  the  data. 
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ADDendix  A:  Tabulated  results 


g-sensitivitv  results  for  lead  frame  design,  20  MHz  fundamental.  0.140”  beam  spacing 


Unit# 

Frequency 

(MHz) 

Q(k) 

R1 

(ohms) 

CO 

(PF) 

Cl 

(fF) 

Hn 

1 

20.061 

80.5 

15.5 

4.0 

6.35 

0.01 

0.01 

0.13 

0.13 

2 

20.006 

91.4 

14.2 

3.8 

6.12 

0.04 

0.23 

0.24 

0.33 

3 

19.996 

79.0 

15.3 

4.0 

6.57 

0.01 

0.27 

0.05 

0.27 

5 

20.088 

82.4 

14.7 

4.0 

6.55 

0.31 

0.06 

0.33 

6 

20.036 

73.8 

17.7 

4.1 

6.08 

0.10 

0.24 

0.09 

0.28 

7 

20.044 

85.4 

14.6 

4.2 

6.39 

0.06 

0.01 

0.14 

mssm 

8 

20.048 

90.3 

13.9 

4.0 

6.34 

0.06 

0.04 

0.05 

0.09 

9 

20.029 

96.0 

13.0 

3.9 

6.39 

0.03 

0.28 

0.05 

0.29 

10 

20.048 

19.7 

67.1 

4.2 

5.99 

0.09 

0.11 

0.12 

0.19 

11 

20.005 

49.1 

26.6 

3.9 

6.09 

0.09 

0.09 

0.12 

0.17 

12 

20.005 

67.2 

17.9 

4.0 

6.61 

0.08 

0.25 

0.09 

13 

20.046 

70.2 

17.2 

4.0 

6.59 

0.08 

0.29 

0.06 

0.31 

MEM 

20.064 

lEaOi 

14.2 

4.0 

6.06 

0.10 

0.04 

0.19 

0.22 

mm 

20.036 

13.9 

4.2 

6.60 

0.15 

0.05 

0.10 

0.19 

mm 

19.995 

95.6 

12.9 

4.0 

6.43 

0.04 

0.22 

0.03 

0.23 

18 

20.054 

80.3 

15.6 

3.8 

6.33 

0.11 

0.14 

0.23 

0.29 

19 

20.013 

69.5 

17.8 

3.9 

6.42 

0.13 

0.11 

0.25 

20 

20.052 

27.6 

44.3 

4.0 

6.49 

0.11 

0.03 

0.10 

MiKEM 

21 

20.023 

72.4 

17.0 

4.0 

6.45 

0.11 

0.17 

0.13 

MiIiEM 

22 

20.032 

83.5 

15.2 

4.0 

6.26 

0.18 

0.16 

0.10 

■«WsM 

23 

16.2 

4.0 

6.54 

0.15 

0.09 

0.22 

MiftW 

24 

20.014 

91.4 

14.0 

3.8 

6.20 

0.07 

0.17 

0.05 

0.19 

25 

20.057 

80.8 

16.1 

3.8 

6.11 

0.17 

0.12 

0.15 

0.26 

26 

20.012 

87.4 

14.2 

6.9 

6.39 

0.04 

0.20 

0.25 

0.32 

27 

20.072 

54.7 

22.7 

6.9 

6.37 

0.04 

0.13 

0.07 

0.15 

28 

20.038 

84.6 

14.8 

4.0 

6.36 

0.06 

0.21 

0.14 

0.26 

30 

20.052 

46.5 

28.0 

3.8 

6.10 

0.05 

0.15 

0.03 

0.16 

31 

20.068 

79.9 

16.0 

3.9 

6.20 

0.05 

0.01 

0.08 

0.09 

32 

20.075 

32.1 

39.3 

3.8 

6.28 

0.08 

0.38 

0.25 

0.46 

33 

20.024 

68.3 

18.2 

3.9 

6.40 

0.05 

0.15 

0.04 

0.16 

34 

20.045 

88.7 

13.8 

4.0 

6.50 

0.08 

0.29 

0.34 

35 

20.012 

63.1 

19.3 

3.9 

6.53 

0.04 

0.16 

0.12 

0.20 

37 

20.020 

81.3 

15.4 

3.8 

6.34 

0.09 

0.12 

0.09 

0.17 

38 

20.071 

58.4 

22.4 

3.9 

6.06 

0.03 

0.10 

0.27 

39 

20.037 

85.0 

14.5 

3.9 

6.45 

0.13 

0.34 

0.18 

0.41 

40 

20.010 

98.9 

12.9 

4.0 

6.25 

0.02 

0.25 

0.05 

MiMsMI 

41 

20.050 

93.8 

13.1 

4.0 

6.48 

0.19 

0.15 

0.07 

Km 

42 

20.065 

59.2 

20.9 

4.1 

6.42 

0.16 

0.15 

0.10 

0.24 

43 

20.033 

60.7 

21.3 

3.9 

6.15 

0.05 

0.03 

0.22 

0.23 

44 

20.048 

13.4 

4.0 

0.08 

0.06 

0.05 

0.11 

45 

20.014 

83.5 

14.6 

7.0 

6.53 

0.15 

0.13 

0.23 

46 

20.052 

73.2 

17.3 

3.9 

6.27 

0.07 

0.04 

0.08 

0.11 

48 

20.031 

82.5 

15.9 

49 

20,039 

106.3 

11.9 

50 

20,078 

69.1 

18.5 

51 

20.055 

40.3 

31.5 

52 

20.045 

82.5 

16.0 

54 

20.047 

74.1 

16.5 

55 

20.049 

92.5 

14.2 

56 

20.019 

85.8 

15.1 

57 

20.043 

59.9 

21.4 

58 

20.072 

36.4 

34.9 

60 

20.063 

83.3 

14.6 

62 

20.024 

51.4 

25.0 

63 

19.987 

55.1 

22.4 

64 

20.057 

87.5 

14.3 

65 

20.035 

69.4 

18.2 

66 

20.024 

77.9 

15.7 

0.09 

0.21 

0.03 

0.13 

0.03 

0.08 

0.11 

0.02 

0.07 

0.10 

0.07 

0.15 

0.09 

0.01 

0.18 

MiVlM 

0.13 

■silH 

0.17 

0.09 

0.16 

0.14 

0.05 

0.09 

0.06 

0.12 

0.08 

0.15 

0.06 

0.18 

0.20 

0.14 

-sensitivity  results  for  lead  frame  design,  60  MHz  3^'*  overtone,  0.140”  beam  soacin 


Unit# 

Frequency 

(MHz) 

Q(k) 

R1 

(ohms) 

102 

59.981 

67.7 

64.8 

103 

59.956 

66.7 

64.4 

104 

59.875 

56.1 

78.1 

105 

59.852 

65.6 

62.3 

106 

59.872 

39.4 

127.4 

107 

59.940 

57.5 

70.6 

108 

59.984 

49.2 

90.2 

109 

59.961 

38.6 

117.1 

110 

59.881 

38.9 

109.7 

111 

59.980 

48.9 

102.8 

112 

59.849 

38.4 

104.1 

113 

59.991 

52.1 

90.9 

114 

59.813 

38.2 

112.1 

115 

59.913 

47.1 

87.6 

116 

59.747 

50.8 

84.9 

117 

59.919 

59.9 

70.1 

118 

59.850 

51.5 

88.6 

120 

59.884 

61.6 

66.6 

121 

59.923 

57.1 

75.3 

122 

59.959 

54.8 

104.5 

123 

59.864 

60.6 

68.8 

124 

59.896 
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Macom  HH-106 

Figure  1.)  Schematic  diagram  of  measurement  system. 


Number  of  Units 


0.0  6.0  10.0  16.0  20.0  26.0  30.0  36.0  40.0  46.0  60.0  66.0  60.0 

Gamma  (10'‘(-10)/q) 

Figure  4.)  Distribution  of  results  -  Z’  elongated,  double-ended  mount. 
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Figure  5.)  Distribution  of  results  -  X’  elongated,  double-ended  mount 
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Figure  6.)  Distribution  of  results  -  Z’  elongated,  single-end  mount 
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Overhang  Distance  (mm) 

Figure  ll.b.)  Variation  of  with  overhang  length  for  various  values  of  support  beam  spacing 
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Figure  13.a.)  Photochemical-Etched  lead  frame. 


Figure  13.b.)  Lead  frame  mounted  Rectangular  AT  Cut  crystal. 


Figure  14.a.)  Results  for  various  lead  frame  mounted  devices. 


Support  Bar  Spacing  (mils) 

Figure  14.b.)  Variation  of  F  magnitude  with  support  beam  spacing,  lead  frame  mount. 
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Figure  16.)  Comparison  of  F  magnitude  for  fundamental  and  3rd  O.T.,  lead  frame  mount. 
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g-senstttvKy  (10-9/g) 


Fig  17a.  Distribution  of  F  in  each  of  the  three  axes,  lead  frame  design, 
20  MHz  fundamental,  0. 140”  beam  spacing 
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Fig  17b.  Distribution  of  |r|,  lead  frame  design,  20  MHz  fundamental, 
0.140”  beam  spacing 
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Fig  17c.  Distribution  of  F  in  each  of  the  three  axes,  lead  frame  design, 
60  MHz  3rd  overtone,  0.140”  beam  spacing 


g-sensttMty  (10-9/g) 


Fig  17d.  Distribution  of  |r  |,  lead  frame  design,  60  MHz  3rd  overtone, 
0.140”  beam  spacing 


58 


-100  -80  -60  -40  -20  0  20  40  60  80  100 

Blank  Orientation  (deg.) 

Figure  18.)  Variation  of  F  components  with  blank  orientation,  SC  Cut,  double  end  mount 
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Figure  19.)  Slotted  clip  for  mounting  rectangular  crystals. 


Figure  20.)  Rectangular  AT  Cut  Z’-elongated  crystal  mounted  with  slotted  clips, 


130  135  140  145  150  155  160 


Clip  Spacing  (mils) 


Figure  21.)  Variation  of  F  with  support  spacing  for  slotted  clip  mounted  crystals. 


Part  3:  Report  on  the  Development  of  the  Miniature  Low  Power  SOIL  &  DDIL 
OCXOs  for  Low  G-sensitivity  applications. 


This  part  consists  of  slides  and  textual  material  that  was  delivered  and  discussed  as 
part  of  the  final  contract  review  meeting. 
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ARMY  CONTRACT#N66001-97-C-8635 
Low  G-Sensitive  Quartz  Resonators,  and 
Low-Power  Clock  Utilizing  the  Resonators 

SOIL  &  DDIL  OCXO  DEVELOPMENT 
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Background  &  Concept 


•  OCXO's  insulation  can  be  improved  by  evacuating  the  package. 
If  aging  performance  is  adequate,  then  an  open  blank  crystal  can 
be  utilized  directly. 

•  Utilize  proven  cold  weld  techniques  for  hybrid  package. 
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The  overall  intent  of  this  development  is  to  produce  a  new  class  of  OCXO’s 
with  smaller  size,  lower  power  dissipation  and  better  characteristics  in 
vibration  than  existing  technology  can  deliver.  The  overall  concept  was 
developed  after  a  careful  analysis  of  previous  efforts  to  develop  evacuated 
OCXO’s.  The  advantage  of  evacuating  the  package  is  a  dramatic 
improvement  in  thermal  insulation.  If  outgassing  and  contamination  issues 
can  be  resolved,  then  the  resonator  does  not  require  a  separate  package  and  the 
entire  module  becomes  dramatically  simpler  (and  smaller).  Achieving 
adequate  cleanliness  was  a  major  concern  from  the  outset  of  this  project  and 
was  the  reason  why  the  aging  testing  was  started  as  quickly  as  possible. 
Although  a  number  of  packaging  options  were  considered,  including  ceramic, 
it  was  ultimately  determined  that  the  best  chance  for  success  was  to  utilize 
cold  weld  sealing.  Cold  welding  is  the  standard  choice  for  high  precision 
resonators,  but  is  essentially  not  used  for  hybrid  microcircuits.  The  Single  DIL 
package  was  previously  tooled  for  a  different  application,  but  the  DDIL 
package  and  sealing  dies  needed  to  be  designed  and  procurred. 
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Material  Selection  &  Processing 


•  BeO  used  (only  in  substrate  form)  for  the  heat  spreading 
substrate 

•  Conventional  thick-film  materials  used  on  the  substrate 


VECTROIM 

international 
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The  above  flow  chart  shows  a  typical  process  flow  for  these  devices.  BeO  was 
selected  to  minimize  thermal  gradients  across  the  temperature  controlled 
substrate,  and  conventional  aluminum  oxide  was  chosen  for  the  lower  level. 
Standard  thick-film  metallization,  epoxies  and  inks  were  used.  Electrical 
components,  including  inductors,  were  chosen  from  Vectron’s  standard  hybrid 
component  inventory.  Although  these  components  have  been  proven  over 
many  years  in  lower  stability  XO’s  and  VCXO’s,  the  outgassing 
characteristics  needed  to  be  proven  in  this  application. 


66 


DIL  OCXO  Platforms 


•  SDIL  OCXO -10MHz 

•  SDIL  OCXO -100MHz 

•  DDIL  OCXO -10MHz 

•  DDIL  OCXO -100MHz 
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The  original  intent  was  to  have  both  HF  and  VHF  devices  with  two  levels  of 
performance,  thus  the  choice  of  10  &  100  MHz  with  two  package  sizes.  HF 
devices  exhibit  superior  aging,  while  VHF  is  an  advantage  when  generating 
microwave  signals.  Since  phase  noise  degrades  by  201ogN  when  the 
frequency  is  multiplied,  lower  multiplication  ratio’s  result  in  lower  phase  noise 
degradation. 
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This  slide  simply  shows  major  elements  of  the  SDIL  EMXO.  This  platform  is 
the  simplest  with  only  one  substrate,  and  all  elements  are  heated.  The 
oscillator  is  essentially  a  CMOS  gate  type  with  added  varactor  diode  and  LC 
trap  for  the  overtone  select.  The  resonator  is  a  3rd  OT  FC  cut.  The  FC  cut 
was  chosen  for  the  following  reasons  1)  better  temperature  performance  than 
AT  cut,  2)  smaller  blank  diameter  for  a  given  motional  Cl  than  SC  and,  3) 
reduced  b-mode  activity  compared  to  SC  cut.  The  primary  difficulty 
encountered  was  eliminating  coupled  mode  activity  around  turnover. 
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As  can  be  seen  in  the  drawing,  the  substrate  is  heated  and  heat  flows  to  the 
resonator  through  the  mounting  clips.  Three  small  posts  are  used  to  suspend 
the  substrate  off  the  package  floor  and  thermally  isolate  the  heated  OCXO 
substrate  from  the  metal  housing.  The  metal  package  is  sealed  with  the  cold¬ 
welding  technique  under  a  high  vacuum  level  environment.  I/O  interconnects 
are  provided  by  wirebonding  the  thermally  isolated  substrate  directly  to  the 
posts  of  the  package.  As  depicted  in  the  above  slide,  the  overall  dimension  of 
the  SDIL  OCXO  is  0.825”x0.520”x0.260”.  The  net  weight  is  slightly  under 
3gm. 
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SDILOCXO.IOOMHZ 
FINCTIONAL  BLOCK  DIACR.\M 


With  the  100  MHz  SDIL  devices,  it  w;is  desirable  to  sink  the  output  amplifier 
heat  to  the  package  case.  Thus,  a  two  -evel  approach  was  developed  wherein 
the  circuitry  which  is  not  temperature  sensitive  is  located  on  the  lower 
substrate.  This  is  critical  since  heat  rise  due  to  uncontrolled  power  can  cause 
the  oven  to  prematurely  shut  down  at  higher  temperatures.  Given  the  very 
high  thermal  resistance  of  around  500 'CAV,  this  becomes  a  more  critical  factor 
than  in  conventional  OCXO’s. 
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SOIL  OCXO-100MHZ 
SIDE  VIEW 


The  lOOMHz  SDIL  design  uses  the  same  package  as  the  10  MHz  SDIL,  but  an 
additional  substrate  is  added  for  the  reasons  described  on  the  previous  slide. 
The  temperature-controlled  portion  of  the  OCXO  is  housed  on  the  top 
substrate  with  the  BeO  material.  The  Bi-level  interconnection  is  joined  via 
wirebonding.  Although  the  overall  dimension  is  exactly  the  same  as  the 
lOMHz  SDIL,  the  net  weight  is  roughly  2  gm  higher. 
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The  oscillator  for  the  10  MHz  DDIL  is  a  Colpitts  configuration  with  a  bipolar 
transistor.  As  expected,  the  Allan  deviation  is  much  superior  on  this  part  as  is 
the  overall  stability.  In  addition,  the  temperature  is  set  with  a  programmable 
pot  thereby  allowing  fine  adjustments  after  sealing.  In  addition,  the  lOMHz 
DDIL  also  utilizes  the  2-level  substrate  design  approach. 
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DDILOCXO-10MHZ 
SIDE  VIEW 


This  design  also  uses  a  two  level  structure  to  dissipate  unregulated  heat.  The 
package  has  also  been  enlarged  to  accommodate  a  larger  crystal  and  more 
sophisticated  circuitry. 
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The  electrical  design  of  the  100  MHz  DDIL  is  also  a  Colpitts,  but  with  the 
necessary  modifications  for  frequency. 
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SOIL  10MHz  CURRENT  AGING 
(Room  Temperature-uncontrolted) 


Measured  heater  current  over  an  extended  time  period.  Changes  in  current 
indicate  an  internal  atmospheric  change.  The  variation  in  average  current 
drain  is  a  result  of  different  oven  temperature  set  points.  Constant  power 
consumption  over  extended  periods  of  time  is  indicative  of  low  leak  rates. 
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Temperature  Stability  Data 
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In  this  and  the  following  charts,  the  best  and  worst  units  are  shown  as  the 
lower  and  upper  X’s  respectively,  and  the  mean  is  the  small  square  in  the 
middle.  The  three  lines  on  the  larger  box’s  represent  25,  50  and  75th 
percentile  distribution  of  performance.  Most  of  the  variation  on  SDIL  units  is 
caused  by  the  resolution  of  setting  the  substrate  temperature  to  the  crystal 
turnover  point.  Although  we  have  identified  some  techniques  to  improve  the 
accuracy  of  setting  the  temperature,  time  contraints  did  not  allow  us  to 
implement  them. 
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Frequency  stability  versus  temperature  of  Single  DIP  units. 

Nominal  frequency  of  10  MHz.  Stability  ranges  from  +/.05  to  0.25  ppm  over 
the  temperature  range  of  -40  to  +70  degC.  The  variation  in  performance  is  due 
mainly  to  the  differences  in  turnover  temperature  of  the  crystals  used.  For  this 
test,  a  standard  oven  temperature  was  used  for  all  crystals,  resulting  in  non- 
optimized  temperature  characteristics.  Individual  adjustment  of  oven 
temperature  will  result  in  superior  performance.  ' 
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As  was  stated  previously  the  best  and  worst  units  are  shown  as  the  lower  and 
upper  X’s  respectively,  and  the  mean  is  the  small  square  in  the  middle.  The 
three  lines  on  the  larger  box’s  represent  25,  50  and  75th  percentile  distribution 
of  performance.  The  temperature  performance  of  the  SDIL  100  MHz  unit 
generally  is  not  as  good  as  that  of  the  10  MHz  SDIL  imits.  This  is  due  in  large 
part  to  a  greater  circuit  sensitivity  at  1 00  MHz  to  changing  stray  reactances 
and  other  circuit  parameters 
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Frequency  stability  versus  temperature  of  Single  DIP  units.Nominal  frequency 
is  100  MHz.  The  crystals  used  in  these  units  are  5th  overtone  AT  cut.  A 
single  modification  was  made  on  each  unit  to  adjust  oven  temperature  to 
match  the  crystal  parameters.  While  tl  ere  is  an  apparent  improvement,  the 
slope  of  the  frequency  /temperature  cu  we  indicates  that  the  oven  temperature 
is  still  set  low  for  most  imits.  Further  individual  adjustment  of  temperature 
and  servo  gain  would  result  in  better  temperature  performance. 
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As  can  be  seen  from  the  data,  these  Double  DIL  units  exhibit  much  better 
stability  than  Single  DIL  units.  This  is  primarily  attributed  to  having  an  ability 
to  adjust  the  oven  temperature  with  the  digital  potentiometer  after  sealing. 
Units  are  run  over  temperature  several  times  at  different  oven  temperature 
settings.  Based  on  the  frequency  temperature  data,  an  optimal  oven 
temperature  is  selected. 
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This  data  indicates  that  most  of  the  frequency  deviation  is  due  to  second  and 
third  order  temperature  effects.  We  believe  that  this  is  caused  by  circuit 
temperature  variations  and  shifting  thermal  gradients.  Both  effects  can  be 
further  improved  with  the  use  of  SC  crystals,  but  at  the  expense  of  tuning 
range  and  possibly  acceleration  sensitivity.  We  do  plan  on  introducing  SC 
cuts  to  improve  thermal  stability. 
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Data  is  shown  to  +65°C  for  the  100  MHz  DDIL  units.  Crystals  exhibited  a 
lower  than  expected  turnover  point,  so  it  was  necessary  to  reduce  the  heated 
substrate  temperature.  Testing  indicated  that  there  was  a  shift  in  apparent 
turnover  from  static  f-t  measurements  of  as  much  as  8°C.  We  suspect  that  this 
is  due  to  thermal  stress  effects  with  the  AT  cut.  Improvements  were  obtained 
by  putting  a  thermal  shield  over  the  crystal,  but  the  temperature  stability  of  the 
10  MHz  devices  was  still  superior.  It  is  felt  that  in  most  applications,  the 
stability  number  achieved  here  are  generally  adequate  since  aging  will  tend  to 
dominate  at  these  frequencies.  However,  it  will  be  necessary  to  improve  the 
upper  temperature  limitation.  We  feel  that  this  can  be  done  through  a 
combination  of  having  a  higher  crystal  turnover  and  improving  the  thermal 
shielding.  Incorporating  FC,  IT  or  SC  resonators  may  also  improve  the 
stability. 
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Again,  it  can  be  seen  here  that  the  linear  effects  are  nearly  absent,  indicating 
good  resolution  on  setting  the  substrate  temperature.  The  ability  to  adjust 
each  oscillator  individually  in  the  DDIL  configuration  yields  much  better 
temperature  stability.  It  is  also  noted  that  temperature  stability  at  higher 
frequency  is  degraded  dueto  an  increased  sensitivity  to  parametric  and  stray 
circuit  changes. 
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SOIL  tOMHz;  Sumix-  State  CtarciM 


The  “flattening”  of  current  at  high  temperature  on  some  units  is  caused  by 
lower  crystal  turnover.  Units  with  +95°C  crystal  turnover  temperatures 
remained  in  linear  operation  to  +85°C.  We  believe  that  the  variation  in 
thermal  resistance  (reflected  in  slope  here)  is  caused  by  varying  degrees  of 
partial  pressure  in  the  can.  Further  work  is  needed  to  determine  the  root  cause 
of  the  variations. 
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Steady  state  current  versus  ambient  temperature  100  MHz.  Single  DIL 
package.  Possible  variations  in  internal  atmosphere  and  set  temperature 
contribute  to  the  power  variations  from  unit  to  unit.  Higher  average  current  in 
the  SDIL  package  at  100  MHz  compared  to  the  SOIL  10  MHz  xmit  is  primarily 
due  to  the  higher  oven  temperature  in  the  100  MHz  SDIL  units. 
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As  expected,  the  average  oven  current  drain  for  the  DDIL  units  is  greater  than 
that  of  the  SDIL  units.  As  previously  mentioned,  the  “flattening”  slope  at 
higher  temperatures  is  due  to  the  ambient  temperature  approaching  the  internal 
oscillator  temperature,  causing  the  oven  to  shut  down. 


Steady  state  current  versus  ambient  temperature.  100  MHz  DDIL  package. 
The  slope  of  current  drain  versus  ambient  temperature  (thermal  resistance) 
remains  constant  over  almost  all  traces  on  this  plot,  indicating  consistent  low 
leak  rates  and  outgasing  inside  the  units.  Again,  the  absolute  power  drain  is 
greater  for  a  DDIL  package  and  a  higher  frequency  design. 
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This  data  was  taken  at  +25°C,  and  represents  measurements  taken  on  10 
pieces.  Individual  warm-up  curves  on  selected  units  is  shown  on  later  slides  at 
various  temperatures. 
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SDILOCXO@100MHz 


This  data  was  also  taken  at  +25°C,  and  represents  measurements  taken  on 
representative  pieces.  Individual  warm-up  curves  on  selected  units  is  shown 
on  later  slides  at  various  temperatures. 

Consistent  with  the  observations  of  stability,  and  power  drain,  the  100  MHz 
units  have  slightly  worse  performance  than  their  10  MHz  counterparts.  Again, 
this  is  primarily  due  to  increased  circuit  sensitivity  at  higher  frequencies. 
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This  graph  showing  warm-up  of  a  representative  group  of  10  MHz  DDIL 
OCXO’s  highlights  the  difference  in  mass  between  the  SDIL  and  DDIL 
packages.  Despite  the  extra  starting  current  used  in  the  DDIL  the  SDIL 
performance  is  still  significantly  better.  The  DDIL  package  contains  a  more 
sophisticated  oven  controller  and  if  faster  warm-up  were  required,  it  would  be 
possible  to  increase  the  oven  startup  current. 
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DDILOCXO@100MHz 


As  might  be  expected,  the  frequency  warm-up  of  the  100  MHz  DDIL  units  is 
the  slowest  of  the  four  types  tested.  In  addition  to  reasons  listed  on  earlier 
pages,  the  transient  response  of  and  AT  cut  crystal  is  significantly  slower  than 
the  doubly  rotated  family  of  crystals,  and  exhibits  significant  overshoot  when 
warming  up. 
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Note  here  that  warm-up  times  can  be  reduced  somewhat  by  increasing  the 
initial  heater  current.  This  is  set  in  the  design  by  the  value  of  one  resistor,  so 
changes  can  readily  be  accommodated. 


96 


Vibration  Sensitivity  Data 
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SDILOCXO@10MHz 
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Statistical  fit  of  r(Gamma) 
Single  DIL  10  MHz  package 
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SDlLOCXO@100MHz 

Histogram 


Absolute  GarTYna(H2/g) 
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Statistical  fit  of  r(Gamma) 
Single  DIL  100  MHz  package 
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DDILOCXO@10MHz 
Histogram 
Absolute  Gamma 


Gaussian  Fit 
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Statistical  fit  of  r(Gamjm.a) 
Double  DIL  10  MHz  package 
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DCHLCX:XO@100(VHz 
Histogram 
Absolute  Gamma 


Statistical  fit  of  r(Gamma) 
Double  DIL  1 00  MHz  package 
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This  table  and  the  associated  data  suggests  that  with  these  resonator  designs, 
100  MHz  5th  overtone  devices  will  exhibit  superior  characteristics.  Since 
there  will  also  be  less  phase  noise  degradation  due  to  frequency  multiplication, 
this  also  suggests  that  100  MHz  is  a  good  choice  when  generating  microwave 
signals  with  good  spectral  purity  during  vibration.  If  better  overall  stability  is 
required  (i.e  temperature  and  aging),  it  would  be  quite  straightforward  to  phase 
lock  the  100  MHz  oscillator  to  a  10  MHz  reference  oscillator. 
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Phase  Noise  Data 


VECTRON 
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A  few  comments  are  in  order  regarding  phase  noise  noise  performance.  All 
units  were  assembled  in  Vectron’s  hybrid  oscillator  facility  in  Norwalk,  CT. 
This  facility  has  been  set  up  to  manufacture  hybrid  XO’s  and  VCXO’s  using 
standard  chip  and  wire  techniques.  Blanks  are  procured  from  outside  suppliers 
and  shipped  to  Norwalk  in  plated  form.  While  this  is  adequate  for  standard 
applications,  it  very  likely  degrades  short-term  stability  and  close-in  phase 
noise  for  this  application.  Consequently,  it  is  our  belief  that  phase  noise  data 
(particularly  variations  at  offsets  of  less  than  1 00  Hz)  may  be  somewhat 
degraded  from  what  ultimately  can  be  achieved.  It  is  Vectron’s  intent  to  either 
find  a  better  method  of  shipping  plated  blanks,  or  establish  baseplating  locally. 
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In  addition  to  the  shipping  issues  previously  commented  on,  we  believe  that 
the  variations  in  setting  the  substrate  temperature  are  impacting  performance 
somewhat.  As  techniques  are  developed  to  increase  the  accuracy  of  setting  the 
temperature,  we  anticipate  improvements  at  small  offset  frequencies. 
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Phase  Noise(dB<yHz) 
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SOILOCXO@100IVHz 


Measured  Phase  Noise 
100  MHz  Single  DIL  Package 


A  voltage  regulator  has  been  positively  identified  as  the  source  of  noise  at 
offsets  of  1 00  Hz  and  larger.  A  new  voltage  regulator  will  be  installed  in 
future  designs. 
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Allan  Deviation  Data 
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Again,  it  is  our  belief  that  better  methods  of  processing  the  blanks  will 
significantly  improved  Allan  deviation.  In  addition,  enhancements  to  the 
accuracy  of  setting  substrate  temperature  on  SDIL  units  may  also  improve 
performance,  particularly  at  10  MHz. 
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SUMMARY 
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R&D  Target 


Size 

inch3 

0.1 

0.125 

4 

Warm-up  Time 

min 

0.5 

0.75 

0.75 

5 

0.1 

Power@25C 

Watts 

0.10 

0.22 

1.5 

Allan  Deviation@1s 

X10-11 

2 

3.84 

100 

Vibration  Sensitivity 

x10-10/g 

1 

6.58 

4.22 

10 

Parameter 

CO-08OCXO 

DDIL  OCXO 

DDIL  OCXO 

OCXO 

TCXO 

Target 

@10MH2 

©lOOMHz 

Typical 

Typical 

Size 

inch3 

1.1 

0.366 

0.366 

4 

0.5 

Warm-up  Time 

min 

2 

1 

1 

5 

0.1 

Power@25C 

Watts 

0.40 

0.33 

0.328 

1.5 

0.05 

Allan  Deviation@1s 

X10-11  Hz 

0.5 

1.55 

8.16 

0,5 

100 

Vibration  Sensitivity 

x10-10/g 

10 

4.99 

3,29 

10 

10 
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Performance  Summary 
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